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Abstract 

High order splitting schemes with complex timesteps are applied to Kolmogorov 
backward equations stemming from stochastic differential equations in Strato- 
novich form. In the setting of weighted spaces, the necessary analyticity of the 
split semigroups can be easily proved. A numerical example from interest rate 
theory, the CIR2 model, is considered. The numerical results are robust for 
drift-dominated problems, and confirm our theoretical results. 
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1. Introduction 

In mathematical finance, the pricing of derivative contracts can be reduced 
to the calculation of expected values under the risk- neutral measure, see, e.g., 
[3I. [38l| . This can be performed in different manners. The most general approach 
is to use Monte Carlo methods, see I2 7| . 1 181. where the recently introduced 
multilevel Monte Carlo approach from [17| led to fast methods for problems 
with very low smoothness. 

In contrast, we want to obtain approximations of the solution of the Kol- 
mogorov backward PDE by splitting schemes. Such methods were also applied 
successfully in Quasi Monte Carlo simulations in [32|], but all such approaches 
are in the end limited by the accuracy of the integration scheme. Furthermore, it 
is not straightforward to evaluate the stochastic processes at complex timesteps, 
which is necessary if splittings of order higher than 2 are to be used. 

We therefore solve the PDE by finite element methods in space and a high 
order splitting method in time. Such high order splitting methods, overcoming 
the order barrier of 2 commonly encountered for splittings with nonnegative 
times, see were introduced in 0, [23| and make use of the analyticity of the 
split semigroups. To show this analyticity, we make use of function spaces en- 
dowed with weighted supremum norms, originally introduced in [35j for proving 
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the existence of solutions to martingale problems for stochastic partial differ- 
ential equations, and used for the analysis of numerical methods for stochastic 



partial differential equations in [12|, [lfj, |13|, [14| . It turns out that using this 
framework, it is very simple to prove analyticity for semigroups stemming from 
stochastic differential equations in Stratonovich form, hence optimal rates of 
convergence follow. In particular, these results apply to problems on unbounded 
domains with unbounded coefficients vanishing at the (finite) boundaries of the 
domain. Such problems are usually difficult to deal with in Sobolev spaces and 
require the use of weighted Sobolev norms vanishing at the finite boundaries, 
see [H, |i[ for some recent results on the Heston stochastic volatility model and 
a thorough discussion of references on this subject. 

A particularly interesting feature of the considered numerical method is that 
the drift part is completely separated from the diffusion part. This means that 
we can choose suitable numerical schemes for each of these parts separately. In 
particular, if the first order hyperbolic drift part is solved by a method such 
as streamline diffusion finite elements or a discontinuous Galerkin method, we 
can expect that the method is robust for vanishing diffusivity. The need for 
such methods, notably for applications in mathematical finance, was recently 
observed in [30(. Other methods yielding such robustness are streamline diffu- 
sion methods, see and discontinuous Galerkin methods, see, e.g., flS 
We stress that an advantage of our method is that different, optimised solvers 
for the drift part and the diffusion part can be used, providing a more flexible 
scheme. 

The paper is organised as follows. Scction[2]recalls the definitions of weighted 



spaces, specialising the results from [35| and [12|, [13j to the finite-dimensional 
case. Furthermore, we show analyticity for a wide class of Markov semigroups in 
the setting of weighted spaces. Next, Section [3] formulates the splitting scheme 
and contains a convergence result. Finally, we show numerical findings for the 
CIR2 model in Section |4j In particular, we observe robustness of the suggested 
method for drift-dominated problems. 



2. The functional analytic framework 

Let us start off by recalling some basic facts from stochastic analysis, see, 
e.g., [H H, H3 for more details. Fix a stochastic basis (f2, T, P, (J-t)t>o) and 
a e?-dimensional standard Brownian motion (W/)j = i d,t>o on it. Let N G N. 
Consider a stochastic differential equation in Stratonovich form on the closure 
D C of a (not necessarily bounded) Lipschitz domain, 

d 

dX(t,x) = ^2Vj(X(t,x))odWi, X(0,x) = x, (I) 
j=o 

where — t and odW^ — At, Vj : D — > M. N are Lipschitz continuous vector 
fields with appropriate smoothness assumptions discussed below in more detail, 
and x <G D. We suppose that the solution (X(t, x))t>o is well-defined in D, in 



2 



particular that up to some common null set N £ J 7 , ¥(N) = 0, X(t, x) £ D for 
all x £ D and t > 0. This is justified by well-known results on stochastic flows, 
see, e.g., [34|, Section V.7]. For / £ C 2 (D) with suitable growth at infinity, Ito's 
lemma shows that u(t, x) :— ¥.[f(X(T— t, x))} satisfies the backward Kolmogorov 
equation 

d -u(t, x) + Cu{t, x) = 0, t > 0, x £ D, (2a) 



df 



u{T,x) = f(x), x £ D. 



(2b) 



Here, C denotes the "sum of squares" partial differential operator, defined for 
g £ G 2 (D) by 



Cg{x) :=V g{x) + \Y,Vfg{x), 



3 = 1 



with Vg(x) :— V{x) ■ Vg(x) the directional derivative for V: D 
g £ C 1 (I?). We split this operator into 

1 d 

C g(x) := V g(x) and Cig(x) := - \fg(x). 

3=1 



The respective split stochastic differential equations read 



-X (t,x) = V a (X°(t,x)), 
at 



X°(0,x)= x, 

dX\t,x) =^2v j (X 1 (t,x))odWjl, X\0,x)=x, 



and 



(3) 
l N and 

(4) 

(5) 
(6) 



3 = 1 



cf. also [32|. 

We recall the following definitions from j35| and 12 , 13, 13] ■ 

Definition 1. A function ip: D — > (0,oo) is called D-admissible weight function 

if limn a .|i_ >clo ijj{x) — oo and ip is bounded on compact sets, where ||-|| denotes 

xeD 

the Euclidean norm. 



Note that the definition from [14 
pact. While the results in 



and 



I simpli 



plifies in this D is locally corn- 

are stated for real- valued functions 



only, they also hold true for the complex- valued versions of the spaces considered 
here. 

Definition 2. Fix k £ No- For j = 0, . . . , k, let ipj : D — > (0, oo) be D-admissible 



weight functions. The space B^ 3 ^ 3 °' 



k {D) is defined as the closure of C%(D), 



the space of functions f:D—>C such that / is bounded and k times differen- 
tiable with all derivatives up to order k continuous and bounded, with respect 



to the norm 



l(^) J= o k ,k, where 

ll/ll(^)j=0,...,fc,*! 



= ll/lk+El-fe 

3 = 1 



(7) 
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with 

U/IU := sup Vo^ 1 l/MI and |/|^., := sup ^- (a:) ~ 1 1 1 2^" / (a;) 1 1 . (8) 
xeD xeD 

Here, f(x) is seen as a j-linear form (C N y — > C, endowed with the norm 

\\D^f(x)\\:= sup \Dlf(x)(z lr -- ,Zj)\. (9) 
zi,..., 2j -ec lv 

|^l | ,- ■ ■ , || Zj I < 1 

We also write B^(D) := Bjf (£>). 

It is shown in [13, El that i3^ j)j=0 '' '"(D}_C C fe (L>), and that it is a Banach 
space. Furthermore, the results of 3f| [l2|, [lOj yield the following generalised 
Feller property. 

Proposition 3. Let ip be a D-admissible weight function. Suppose that (Pt)t>o 
is a family of bounded linear operators on Bq(D) such that 

• Pq = I , the identity on B^(D), 

• Pt+s = PtPa for all t, s > 0, and 

• lim^o Ptf(x) = f(x) for all x G D and f e B^'(D). 

Then, (Pt)t>o is strongly continuous on B^(D). 

For simplicity, we formulate the following result for polynomially growing 
functions and vector fields with bounded derivatives. Therefore, we fix s G N 
large enough and set tpe(x) '■= (1 + H^H 2 )^ - ^ 2 , f € No and ip(x) := ipo(x). 
Similar results can be obtained for different choices, see, e.g., [l3| for the case of 
exponentially growing functions and bounded vector fields. It is straightforward 
to see that P t f(x) := E[f(X(t,x))], P t °f(x) := E[f(X°(t,x))] and P£f(x) := 
E[f(X 1 (t, x))] define strongly continuous semigroups on B^(D), and that their 
generators are suitable extensions of £, Co and Ci, respectively. Here, we apply 
the usual time inversion to turn the final value problem @ into an initial value 
problem. 

Proposition 4. Let k e N. Assume that for j = 0, . . . , d, Vj € C k+1 (D:R N ) 
with all derivatives of order 1 to k + 1 bounded. Then, P t G^(D) C C^(D) for 
all t > 0, and 

11^/11(^=0,...,* < C{t)\\f\\ m . =0 k for all f G B^= - ' k (D), (10) 

where C: [0,oo) — > (0,oo) is monotonously increasing. In particular, 

P t : B^ j)i=0 '-' h (D) -> B^ ih=0 -- k (D) (11) 
is a bounded operator with norm < C(t), t > 0. 
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Proof. The first part and the norm bound is a consequence of [28j, Theorem 3.3, 
p. 223], and the invariance of B^ °""' k (D) follows by a density argument. 
See also [l3|, Lemma 13] for a related result. □ 

Proposition 5. Let V: D — > K be a continuously differentiable vector field 
with bounded derivative. Suppose that V(x) ■ n(x) < for almost every x £ dD, 
where n: dD — » R w is the (almost everywhere defined) outer unit normal vector 
to D. For t > 0, consider the semiflow Fl^ : D — > D ofV, i.e., 

±F\Y(x)=V(F\Y(x)), F%(x)=x. (12) 

Then, 

PYf{x):=f{F\V{x)), t>0, (13) 

defines a strongly continuous semigroup on B^(D), and its generator is the 
closure of the operator V : B^'^ (D) -» B 4 '{D). IfV{x) ■ n(x) = on dD, 
(P^)t>o can be extended to a strongly continuous group. 

Proof. The Lipschitz continuity of V together with V(x) ■ n(x) < for almost 
every x £ dD yield that Fl t is continuous and well-defined for t > 0, and that 
for some e > there exists C > such that 

ip{FlY(x)) < Cty(a?) and (14) 
i>i{F\Y (x)) < Cipi(x) for all t £ [0,e) and x € D. (15) 

Hence, for arbitrary / £ B^(D), 

sup^r 1 !/^*))! < sup^(a:)-V(Fl t y (a ; ))||/|U,o < C||/|U,o, (16) 

x£D x£D 

whence P t v : B^(D) B^(D) is bounded for t > 0. As the semigroup property 
of (P t y ) tS R is obvious, we can apply Proposition[3]to obtain that (P t v ) t >o indeed 
is a strongly continuous semigroup on B^(D). 

Let us now determine its generator Q v . We shall prove that b[^'^ 1 \d) C 
dom£ y , and that G v f(x) = Vf{x) on this space. Clearly, for / £ C£(P>), 
^| t=0 /(Fl t y (x)) = Vf(x), and V f £ B*(D). Hence, 

PYf-f=fpY{Vf)ds, (17) 

and this implies that / £ domQ v and Q v f = V f . The Lipschitz continuity of 
V implies ||V(a;)|| < C(l + Hxlj 2 ) 1 / 2 with some constant C > 0, and we observe 

\vf(x)\ < c(i + IMI 2 ) 1/2 <Mz)l/k,i < c^(x)\fy ul . (is) 

It follows that V: B^'^iD) -» B^{D) is bounded, and together with the 
closedness of Q v and Propostion [TT] this implies Bi'(D) C AouyQ v and 

g v f = Vf for all / e B^'CD). (19) 
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we see that x h-» P t v f(x) is continuously differentiablc, 

and that for z £ C , 

\DPYf(x)z\ = \Df(F\Y(x))DF\Y(x)z\ < ^^(Ftf (x))\DF% (x)z\. (20) 
Denote the Lipschitz constant of V by L, then by Gronwall's inequality, 

\DF\Y(x)z\<exp(L\t\)\\z\\. (21) 

Hence, fort G [0,e), \\P^f\\ { ^ l)A < C\\f\\ { ^ lhl . As additionally P t v (C£(D)) C 

CJ(D), we obtain that B^^\D) is a core for Q v by [15, Proposition II.1.7]. 

If we assume that V(x) ■ n(x) — for x £ dD, we can apply the above results 
to both V and — V, from which it follows that (P t y ) t6 R is a strongly continuous 
group. The proof is thus complete. □ 

Corollary 6. Let V £ C 3 (D;1$L N ) with bounded derivatives of order 1 to 3 such 
that V{x) ■ n{x) — for almost every x £ dD. Then, 

P? V2 f(x) := ^/(Fl^(x))^Lexp(-^)d y (22) 

defines an analytic semigroup of angle ir/2 on B^(D). Its generator is the 
closure of the operator \V 2 : B^^ 1 ^ (D) -> B*(D). 



Proof. This is a consequence of Proposition [5] and 15, Corollary II. 4. 9]. That 

i 4 



B^^ 1, ^ 2 \d) is a core is a consequence of Proposition |H □ 



Corollary 7. Suppose that the vector field Vo is continuously differ entiable 
with bounded derivative and that the vector fields Vj, j = l,...,d are three 
times continuously differ entiable with bounded derivatives of order 1 to 3 on the 
Lipschitz domain D C M. , and satisfy 

Vo(x) ■ n(x) < and Vj(x) ■ n(x) = for j — l,...,d and x £ dD, (23) 

where n(x) denotes the (almost everywhere defined) outer unit normal vector 
on D in x £ dD. Define C^g{x) := Vog(x) and C\g{x) := ^ Sj=i V?9( x ) f or 

g £ B^' 3 ^ 3=0 ' 1 ' 2 (D). Then, the closure of Co generates the strongly continuous 
semigroup (P t °)f>o, and the closure of C\ generates the analytical semigroup of 
angle tt/2 (P£) t >a on B^{D). 

Proof. This follows from a combination of Proposition [5l Corollary [6j and 
Proposition [T2l □ 

Remark 8. We only obtained properties of the split semigroups [P[ )t>o and 
(Pt)t>o in Corollary [7] While perturbation theory arguments should allow 
us to derive properties of the semigroup (Pt)t>o from these results by applying 
perturbation theory, see, e.g., [15|, Theorem III. 2. 10], this question is not central 
for the topic of this paper and is therefore left to further research. 



G 



3. High order splitting schemes 

Recently, it was observed in 2^, g| that complex timesteps can be used to 



obtain splitting schemes with order higher than 2 for parabolic problems. In 
the given setting, we have the following result. 

Theorem 9. Let (<Xk)k=i,...,s C [0,oo), (Pk)k=i,...,a C {z £ C: Rez > 0} be 
the coefficients of a splitting scheme of formal order p. Assume the setting 
of Corollary [?] and that additionally, Vj for j = 0, . . . , d are 2p + 3 times 
continuously differentiable with all derivatives of order 1 to 2p + 3 bounded. For 
f € B'^(D), set 

Q(At)f ~ P^AtP^At ■ ■ ■ pIaiPIaJ- (24) 

Then, for f € ^(p+i) ' "' 2(p+1) (D) and T > 0, there exists a constant Ctj > 
independent of n £ N such that 

\\Q n ( T,n)f-PTfh<C T jn-v. (25) 
Proof. By Proposition H b(^-'*»c»+u) 

(D) is invariant with respect to the 
semigroup (Pt)t>o- Using the results from Section El th is is hence a straight- 



forward consequence of [22|, Theorem 2.3], see also [23J, Theorem 3.1] and 



Theorem 3.1]. □ 

Remark 10. A fourth order splitting as required by the above theorem was 
constructed in 0, equation (5.1)]. They have s — 5 and 

a% = 0, et 2 = a 3 = a 4 = a 5 = -, (26) 

This scheme requires the evaluation of the semigroups P^ k A* anc ^ Pp k At t° fourth 
order. The strategy of composition methods as applied in [6j will fail if Cq does 
not generate a group, which is a consequence of As Vo usually does not 
satisfy Vq(x) ■ n(x) — for x £ dD in problems from mathematical finance 
with mean reversion because Vo points towards the mean and hence Cq does 
not generate a group, we do not analyse this question further here. 



4. Numerical example: The CIR2 model 

The CIR2 (Longstaff-Schwartz) model, see [HQ, is a popular model for 
interest rates. It supposes that the short rate r t is given by 

n = x t + y t (28) 

with xt and yt being square-root diffusions (CIR processes) under the risk- 
neutral measure Q, 

dx t = 9 x (p x - x t )dt + o- x -s/xl6Wf, (29a) 
dy t = 6 y (fi y - y t )dt + o-yy/yldW? (29b) 
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where the stochastic basis is given by (fi, J 7 , Q, (Ft)t>o), d x , fi x , a x , 9 y , fi y and 
a y are positive constants, and Wf and Wf are independent Brownian motions. 
From (r t ) t >o, prices of default-free zero coupon bonds can be derived via 



B(S,T) :=E[exp(- / r t dt)\T s 
Js 



(30) 
and 



where E[-|J r 5] denotes conditional expectation with respect to Fs under 
< S < T. 

For simplicity, we consider the approximation of the bond price itself, for 
which analytic formulas are available. This allows us to determine the precise 
error of our numerical results. As the bond price has smooth dependence on 
the initial values xq and yo > we do not need to use non- uniform timestcps in the 
approximation. In the approximation of non-smooth payoffs, however, this is 
necessary; for a construction of non- uniform time grids, see, e.g., 

After the usual time inversion, the PDE formulation of the problem reads 



d_ 

dt 



u{t,x,y) = Cu(t,x,y), u(0,x,y) = l, 



(31) 



with u(T, x, y) = B(0, T) if we fix xq — x and yo = y, where 
Cu(t,x,y) = 



. J . 1 2 d 2 . 



d Id 2 

ii^y - x , v) + o°"»2/-hz9 x > y)~ ( x + y) u (t, x , v)- 



dy 

We split this equation according to 

1 o 



C u {t,x,y) 



Ciui(t,x,y) 



X i^f^X % 



y) 



K 2 



1 2 d 2 , 



d_ 

dx 
d_ 
dy 
1 



dy 2 

u (t,x,y) 
uo(t,x,y) and 



(32) 



(33) 



1 



d 



-a-- Ul (t,x,y) 
1 



'dy 



(34) 



which corresponds to a split into drift £ and diffusion C\ after transformation 
into Stratonovich form. In order to ensure that x = and y — are outflow 
boundaries for the hyperbolic problem ^u(t,x,y) — Cou(t,x,y), we assume 
^Sx[i-x > &x and 4:9yfj, y > a 2 , which is weaker than the Feller condition, see, 
e.g., 0. While the given operator does not satisfy the assumptions of Section [2] 
to prove that C\ generates an analytic semigroup, this can be proved directly. 
Proposition[T3lshows this result for the ID case, and the proof generalises to the 
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considered equation. The necessary smoothness of the exact solution u(t, x, y) of 
equation (|3ip can also not be obtained directly from Proposition |4l but follows 
from [1, Proposition 43]. 



I le-05 



I le-05 




liilk'ML'p:. 

(b) e = .125 



Figure 1: Errors in approximation of bond prices of the CIR2 model 

In the computations, we fix the parameters 

6 X = 15.5, (i x = .025, a x = .2s, 6 y = 20.5, \i y = .025, a y = .3e, (35) 

where we consider s = 1. and e = .125. In both cases, the problem is drift- 
dominated. As we want to focus on the time discretisation error, we minimise the 
impact of the domain truncation by cutting of far beyond any reasonable values 
for the state variables, at x = 16 and y = 16, and of the spatial discretisation 
by using higher order continuous finite elements in space, with 17 grid points 
per direction and a polynomial degree of 4. This leads to a total of 4225 spatial 
degrees of freedom. The operator splitting uses the scheme from [6J, equation 
(5.1)], see Remark [TU] for the coefficients. For the split equations, we employ 
a streamline diffusion FEM [25] using a polynomial degree of 4 for the first 
order hyperbolic equation involving £q , and a resolvent Krylov subspace method 
using a 30-dimensional Krylov subspace for the approximation of the matrix 
exponential 21, ]j| of the second order parabolic equation involving C\. 
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Figure 2: Effects of domain truncation 

Figure [1] shows the results of a numerical computation. The plotted error 
is the difference between B(0,T) and the numerical approximation i3 app (0,T), 
measured in the £>^-norm with ip(x,y) = (1 + x 2 + y 2 ) 3 . We clearly observe 
4th order convergence. If we restrict x + y < 1, which should include all values 
of practical interest, we see that ip(x,y) < 64. Hence, in this region, the form 
of the error norm ensures that the pointwise approximation error is less than 
10~ 3 using only 8 timesteps. Furthermore, the error is virtually independent 
of the value of e, which means that our approximation is robust for the case of 
vanishing diffusion. 

In order to understand the effect of the space truncation, we plot in Figure [2] 
the errors for different choices of the cutoff point, where the spatial mesh width 
is kept constant and 8 timesteps are used. We observe that if the domain size 
is increased, the error stays bounded. 

Appendix A. Functional-analytic results 

In this section, we collect a result on closed operators and one on analytic 
semigroups that are easy to prove, but that were not able to find in standard 
literature. 

Proposition 11. Let (X, \\-\\x) and (Y, ||-||y) be Banach spaces such that Y is 
continuously embedded in X, i.e., Y C X and \\y\\x < C||y||y for all y € Y 
with some constant C > 0, and Z be dense in Y (with respect to the norm 
ofY). Given a closed operator A: domA C X — > X with Z C domj4 such 
that HA^Hx < ^Iklly f or all z ^ Z with K > 0, we have Y C domA and 
\\Ay\\x < K\\y\\ Y for ally £ Y . 

Proof. By continuity, A\z can be extended to a continuous linear operator 
B.Y^X with ||By||x < K\\y\\ Y and Bz = Az for all z € Z. Given y e Y, 
let (z„) n6 N be such that limn^ooll?/ — z n \\Y = 0. Then, by continuity of B, 
lim„^ 00 ||.Bi/ — Bz n \\x — 0. As Bz n — Az n for all n 6 N, we observe z n — > y 
and Az„ — > By in the norm of X. Hence, the closedness of A implies y € domA 
and Ay = By, which yields the claim. □ 




10 



Proposition 12. Fix a complex Banach space X. Given 6 € (0,7r/2], set 
"Eg '■= {z£C: |argz| <S}. Let (A,domA), (B,domB) be the generators of 
analytic semigroups (S^) Z £-£ 6 uo an d (S z ) z es s uo on A satisfying 

\\Sf\\ < exp(K\z\) and ||Sf || < exp(K\z\) for z e E 5 U (A.l) 

for some K > 0, and assume that the closure (C, domC) of (A + B, dom A (1 
dom_B) generates a strongly continuous semigroup (S^)t>o- 
Then, (Sf)t>o extends to an analytic semigroup on £5 U 0. 

Proof. We assume without loss of generality that the semigroups generated by 
A and B are bounded on £5; otherwise, translate the operators appropriately 



first. By 15|, Theorem II. 4. 6], we have to prove that for some 9 g (0, 7r/2), e ±l6 C 
generate strongly continuous semigroups. The Hille-Yosida theorem implies 
existence of wo € K such that A — C is surjective for Re A > cjo- As Re e ±t8 = 
cos(9) > 0, we can find [i > such that Re^e ±l8 > w , whence fi — e ±l9 C — 
e ±t6 (e Tie ji — C) is invertible. Fixing 9 e (0, 6), we can apply the Trotter product 
formula [lja, Corollary III. 5. 8] to e ±z9 A and e ±l0 B to prove that e ±t9 C generate 
strongly continuous semigroups. Repeating this argument with A, B and C 
replaced by e la A, e la B and e la C, where a e (—6,6), we obtain the result. □ 

Appendix B. Analyticity of the CIR semigroup in the B^ setting 

The following result proves the analyticity of the CIR semigroup on a B^- 
space directly, as Corollary [5] does not apply due to the lacking Lipschitz conti- 
nuity of the square root. 

Proposition 13. Set ip(x) := (1 + x 2 ) 8 / 2 with s e N. For f e B^([0,oo)), 
define 

P t f{x):= Jj{{^+aVty) 2 )^=eM^)^y, *€[0,ao), (B.l) 

Then, P t f defines an analytic semigroup on B^'([0, 00)) and solves the PDE 
d d 2 d 

—u(t,x)=2a 2 x——^u(t,x)+a 2 —u(t,x), t>0andx>0, (B.2a) 
ot ox z ox 

u(0,x) = f(x), x>0. (B.2b) 

Proof. It is easy to see that (Pt)t>o defines a generalised Feller semigroup on 
B^([0, 00)), hence it is strongly continuous by Proposition [3] Integration by 
parts proves 

|p t /(x) = jf f(W* + aV-tyf) ^ ±^ ay-}= exp(-^)d y 

= hl y|-/((^ + ^2/) 2 )^=exp(-^)dy 
2i J R oy V27T 2 

, ); / (y*-l)f((^ + aViy) 2 )-l=e X p(-^)dy, 
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and as the right hand side defines a bounded operator on fi^([0, oo)) with norm 
of the order 0(i _1 ), (Pt)t>o is an analytic semigroup. 
To show that P t f satisfies the given PDE, we calculate 



2^ 



A= / y.f((^ + ^y) 2 )^=exp(-^)dy 

T^tX JR V27T ^ 

and, similarly, 

= i (2^// /(( ^ + ^ )2) vfe cxp( 4) d ^ 



t/ R y l /((V5+ff ^ )2) ^ cxp( 4^ 

TP 1 vf{{ ^ + <7 " / * y)2) vfc exp( ~l" )dy 



4(7 
] 

4cr 2 ta; 

4(7 

] 

4cr 2 t:r 



Hence, 



2 ° 2x & + ° 2 l) p ^ - k jy - i)/((^ + ^) 2 )^ox P (4)d y 

which proves the claim. □ 
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